Effects of fluctuations and Coulomb interaction on the transition temperature of 

granular superconductors. 
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We investigate the suppression of superconducting transition temperature in granular metallic 
systems due to (i) fluctuations of the order parameter (bosonic mechanism) and (ii) Coulomb re- 
pulsion (fermionic mechanism) assuming large tunneling conductance between the grains gr 2S> 1 . 
We find the correction to the superconducting transition temperature for 3d granular samples and 
films. We demonstrate that if the critical temperature T c > grS, where 5 is the mean level spacing 
in a single grain the bosonic mechanism is the dominant mechanism of the superconductivity sup- 
pression, while for critical temperatures T c < grS the suppression of superconductivity is due to the 
fermionic mechanism. 

PACS numbers: 74.81. Bd, 74.78.Na, 73.40.Gk 



I. INTRODUCTION 

Being an experimentally accessible electronic system 
with the tunable parameters granular supercon- 
ductors, offer a unique testing ground for studying com- 
bined effects of disorder, Coulomb interactions and super- 
conducting fluctuations that govern the physics of dis- 
ordered superconductors and are central to mesoscopic 
physics. One of the fundamental questions long calling 
for investigation is the problem of suppression of the su- 
perconducting critical temperature, T c , in granular su- 
perconductors and the role played in this suppression by 
the Coulomb repulsion and superconducting fluctuations. 
In this paper we present a quantitative theory of the sup- 
pression of T c in granular samples. 

The customary belief was that - according to the An- 
derson theorem^- disorder leaves critical temperature of a 
superconductor intact. However this result holds only in 
the mean field BCS approximation, and in all the cases 
where the extension beyond the BCS approximation is 
required, one can expect a noticeable suppression of the 
critical temperature. 

The main mechanisms of the superconductivity sup- 
pression are Coulomb repulsion and superconducting 
fluctuations. For example, disorder shifts significantly 
the superconducting transition temperature in the 2d 
thin films2i2ii2iiiii£. The physical reason for the sup- 
pression of the critical temperature is that in thin films 
the interaction amplitude in the superconducting chan- 
nel decreases due to peculiar disorder-induced interfer- 
ence effects which enhance the effective Coulomb inter- 
action. On the technical side, in order to evaluate the 
effect of disorder, one should sum a certain class of di- 
agrams that include, in particular, cooperons and diffu- 
sons. In the subsequent discussion we will be referring to 
this mechanism of the superconductivity suppression as 
to the fermionic mechanism. 

The superconducting transition temperature can also 
be reduced by the fluctuations of the order parameter, 



the effect being especially strong in low dimensions. The 
corresponding mechanism of the superconductivity sup- 
pression is called the bosonic mechanism. In particular, 
the bosonic mechanism can lead to the appearance of 
the insulating state at zero temperature. The physics 
of this state can be most easily understood in the case 
of a granular sample with weak intergranular coupling: 
the Cooper pair can be localized on a single grain if the 
charging energy is larger than the Josephson energy cor- 
responding to the intergranular coupling—. Later it was 
shown 14 , that a similar mechanism of Cooper pair local- 
ization appears even in the case of the homogenously dis- 
ordered films and the superconductor to insulator tran- 
sition was predicted to occur at zero temperature. 

In this paper we study the corrections to the supercon- 
ducting transition temperature in granular metals per- 
turbativelly. While this approach is restricted and can- 
not be used for study of non-perturbative effects such 
as the superconductor to insulator transition, it is useful 
in a sense that both relevant mechanisms of the critical 
temperature suppression can be studied systematically 
within the same framework. The power of the pertur- 
bative calculation in the study of granular metals was 
demonstrated in that it revealed an important energy 
scale r = g T 6, which was missed for example by the ef- 
fective phase functional formalism, where g T is the tun- 
neling conductance between the grains and S is the mean 
energy level spacing for a single grain, appearing in gran- 
ular materials. The presence of this energy scale which 
has a simple physical interpretation of an inverse aver- 
age time that an electron spends in a single grain before 
tunneling to one of the neighboring grains 15 , brings into 
play new behaviors that are absent in homogeneous me- 
dia. In particular, the two different transport regimes at 
high, T > r, and low, T < V, temperatures appear^ In 
the high temperature regime the correction to the con- 
ductivity due to the Coulomb interaction depends log- 
arithmically on temperature in all dimensions, while at 
low temperatures the interaction correction to conduc- 
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tivity has the Altshuler-Aronov formii and, thus is very 
sensitive to the dimensionality of the sample. 

In a view of these findings one may expect that the 
correction to the superconducting transition temperature 
can also be different depending on whether the temper- 
ture is larger or smaller than the energy scale T. 

In the present paper we analyze the mechanisms of 
the suppression of superconductivity in both temperature 
regimes. We find that the fermionic mechanism is tem- 
perature dependent and that its contribution is strongly 
reduced in the region T > T. In this regime the bosonic 
mechanism of the T c suppression becomes dominant. In 
the low temperature regime, T < T, the correction to the 
critical temperature is similar to that obtained for homo- 
geneously disordered metals. In this regime the fermionic 
mechanism plays the major role as long as the intergran- 
ular tunneling conductance is large. 

The paper is organized as follows: in Sec. ^ we sum- 
marize the results for the suppression of superconducting 
transition temperature of granular metals. In Sec. lIIII wc 
compare our results for the suppression of superconduc- 
tivity in granular metals with known results for homo- 
geneously disordered systems. In Sec. |W] we introduce 
the model; the effect of fluctuations and Coulomb inter- 
action on the superconducting transition temperature is 
then discussed in Sec. El The mathematical details are 
relegated to the Appendixes. 



II. SUMMARY OF THE RESULTS 

It is convenient to discriminate corrections due to 
bosonic and fermionic mechanisms and write the result 
for the suppression AT C of the superconductor transition 
temperature in a form 



AT C / AT C 



T, 



( AT C 



(la) 



where the two terms in the right hand side correspond to 
the bosonic and fermionic mechanisms, respectively. The 
critical temperature T c in Eq. (|la(l is the BCS critical 
temperature. 

We find that at high temperatures, T > T, the 
fermionic correction to the superconducting transition 
temperature does not depend on the dimensionality of 
the sample 



AT C 



-ci 



/ 



d = 2,3. 



(lb) 



where ci — 7£(3)/27r 2 — (ln2)/4 is the numerical coef- 
ficient and d is the dimensionality of the array of the 
grains. 

In the low temperature regime, T < T, the fermionic 
mechanism correction to the superconducting transition 
temperature depends on the dimensionality of the sample 



and is given by 



^-ln 2 ^_ 
l i n 3 r_ 

24 7T 2 g T 111 T c ' 



d = 3 
d = 2 



(lc) 



where A = g T a 3 J d 3 q / (2ir) 3 e ^ 1 m 0.253 is the dimen- 
sionless constant, a is the size of a single grain and 



£ q = 2 9t^(1 - cosqa) 



(Id) 



with {a} being the lattice vectors (we consider a periodic 
cubic lattice of grains). Note that in the low tempera- 
ture regime T < T the correction to the critical temper- 
ature in the dimensionality d = 2 coincides with that 
obtained for homogeneously disordered superconducting 
films upon the substitution T — > r . 

On the contrary, the correction to the transition tem- 
perature due to the bosonic mechanism in Eq. I|la|) re- 
mains the same in both regimes and is given by 



AT C 
Tc 



top J_ 
rem l n s^l 
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d= 3 
d = 2 



(le) 



where C( x ) is the zeta- function and the dimensionless 
constant A was defined below Eq. ((Tc|l . Note that the 
energy scale T does not appear in this bosonic part 
of the suppression of superconducting temperature in 
Eq. Ulejl. This stems from the fact that the character- 
istic length scale for the bosonic mechanism is the co- 
herence length ^ which is much larger than the size of a 
single grain. The result for the two dimensional case in 
Eq. I|lc|) is written with a logarithmic accuracy, assuming 
that \n{gl8/T c ) > 1. 

The above expression for the correction to the tran- 
sition temperature due to the bosonic mechanism was 
obtained in the lowest order in the propagator of super- 
conducting fluctuations and holds therefore as long as the 
value for the critical temperature shift given by Eq. ijle|) 
is larger than the Ginzburg region (AT) G 



(AT) C 



i 

<JT 



d = 3, 
d = 2. 



(2) 



Comparing the correction to the transition temperature 
T c given by Eq. fle| with the width of the Ginzburg 
region, Eq. J2J), one concludes that for id granular metals 
the perturbative result ((Te)l holds if 



T c < g 2 T S. 



(3) 



In two dimensions the correction to the transition tem- 
perature in Eq. (|lef> is only logarithmically larger than 
(AT)g in Eq. (J2J. With the logarithmic accuracy we 
note that the two dimensional result ifTej l holds in the 
same temperature interval @ as for the three dimen- 
sional samples. 
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Note that inside the Ginzburg region the higher or- 
der fluctuation corrections become important. More- 
over, the non perturbative contributions, in particular, 
the contributions from superconducting vortices should 
be taken into account as well. These effects destroy the 
superconducting long range order and lead to Berezinskii- 
Kosterlitz-Thouless transition in 2d systems. 

To summarize our results, we find that the correction 
to the superconducting transition temperature of gran- 
ular metals comes from two different mechanisms; the 
dominant mechanism depends on temperature range. 

(i) In the low temperature regime, T < T, the fermionic 
mechanism is the main mechanism of the suppression of 
T c and the correction to the transition temperature is 
given by Eq. (|lc|) . 

(ii) In the high temperature regime, T > F, the domi- 
nant mechanism is bosonic. At moderate temperatures, 
T < g^S, the correction to the transition temperature 
is perturbative and is given by Eq. (TeJ, while at higher 
temperatures T > g~,6 the superconducting transition 
temperature must be determined by considering the crit- 
ical fluctuations in the effective Ginzburg-Landau func- 
tional. 



III. COMPARISON OF THE SUPPRESSION OF 
SUPERCONDUCTIVITY IN GRANULAR AND 
HOMOGENEOUSLY DISORDERED SYSTEMS 

In this section we compare our results for the sup- 
pression of superconductivity with the known results 
obtained for homogenously disordered superconductors. 
We begin our discussion with homogeneously disordered 
samples briefly reminding what is known about suppres- 
sion of superconductivity in this case. 

Both mechanisms of the suppression of superconduc- 
tivity in homogeneously disordered films were discussed 
in several publications&^iiSiiiiii In particular, for films 
with thickness d such that I <C d <C £ where I is the elec- 
tron mean free path and £ is the coherence length it was 
shown that the result for the suppression of supercon- 
ducting critical temperature can be written in analogous 
form with Eq. (JTaJ, AT C /T C = (AT C /T c ) f + (AT C /T c ) b , 
where^ 



ln 3 [l/(rT c )], 



and 



( AT C 



c / b 



7C(3) 
2^g 



H9/(rT c )}. 



(4a) 



(4b) 



Here g ^> I is the film conductance (per one spin compo- 
nent) and r is the elastic electron mean free time. One 
can see from Eqs. Q that in the regime of large con- 
ductance within the logarithmic accuracy the fermionic 
mechanism is the dominant one. At the same time, if 
the conductance is not too large both corrections be- 
come of the order of one and the bosonic mechanism 



becomes very important as well. In this regime the sup- 
pression of superconductivity should be considered non- 
perturbativelly as in Ref. [Ij for the bosonic mechanism 
and in Ref. 40 for the fermionic mechanism. 

In granular superconductors situation is different due 
to appearance of the energy scale T — g T 5. As one can 
see from Eqs. Q both mechanisms of the suppression 
of superconductivity are important. In the limit of high 
temperatures T > T the interference effects in granu- 
lar metals are suppressed and that is why the fermionic 
mechanism is strongly reduced, Eq. (|l b|) . The shift of the 
superconducting critical temperature in this region is de- 
fined by the bosonic mechanism and has a classic origin. 
In the low temperature limit T < V quantum interfer- 
ence effects become important therefore the suppression 
of superconductivity is defined by the fermionic mecha- 
nism. The fact that in the low temperature regime the 
correction to the superconducting transition temperature 
for a granular samples can be obtained from the corre- 
sponding result for the homogenously disordered samples 
via the substitution of the effective diffusion coefficient 
D = gxSa 2 , suggests that the universal low temperature 
description proposed in Ref. can be generalized to in- 
clude the superconducting channel. 



IV. THE MODEL 

Now we turn to the quantitative description of our 
model and derivation of Eqs. Ijlal [Tel ITe|) . We consider 
a d— dimensional array of superconducting grains in the 
metallic state. The motion of electrons inside the grains 
is diffusive and they can tunnel between grains. We as- 
sume that if the Coulomb interaction were absent, the 
sample would have been a good metal at T > T c . 

The Hamiltonian of the system of the coupled super- 
conducting grains is: 



H — H o + H c + H t 



(5a) 



The term Hq in Eq. (|5a|) describes isolated disordered 
grains with an electron-phonon interaction 



Ho — ^ £i,ka\k a i,k — A ^ a\ k a\ _ k ai t -k>cn t k' + H imp , 

(5b) 

where i labels the grains, k = (k, f ), — k = (— k, I); A > 
is the interaction constant; ajt(<^j,fc) are the creation 
(annihilation) operators for an electron in the state k of 
the i-th grain, and Hi mp represents the elastic interaction 
of the electrons with impurities. The term H c in Eq. I|5a() 
describes the Coulomb repulsion both inside and between 
the grains and is given by 



(5c) 
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where Cy is the capacitance matrix and n, is the operator 
of the electron number in the i-th grain. 

Eq. lfKc)l describes the long range part of the Coulomb 
interaction, which is simply the charging energy of the 
grains. The last term in the right hand side of Eq. I|5a|) 
is the tunneling Hamiltonian 



Ht — X! ti 3 a i,p a 3><l> 



(5d) 



where ijj is the tunneling matrix element corresponding 
to the points of contact of i-th and j— th grains and p, q 
stand for the states in the grains. 

In the following section we will study effects of fluctu- 
ations on the superconducting transition temperature of 
granular metals based on the model defined by Eqs. J5j. 



V. EFFECTS OF FLUCTUATIONS AND 
COULOMB INTERACTION ON TRANSITION 
TEMPERATURE 

The superconducting transition temperature can be 
found by considering corrections to the anomalous Green 
function F due to fluctuations of the order parameter 
and Coulomb interaction in the presence of infinitesimal 
source of pairs A? 8 - Without account of fluctuations and 
interaction effects, the anomalous Green function F is 
given by the expression-^ 



F^e n ) = A/(sl+e 



(6) 



where £ = p 2 /2m — fi, and e n = 2irT(n + 1/2) is the 
fermionic Matsubara frequency. The suppression of the 
transition temperature T c is determined by the correction 
to the function F(£,e n ) 



(7) 



where the function SF(^,e n ) represents the leading or- 
der corrections to the anomalous Green function F(£, e n ) 
due to pair density fluctuations and Coulomb interac- 
tion. The function SF(^,e n ) can be found by means of 
two different methods which lead to identical results: (i) 
solving the Usadel equation with the help of perturba- 
tion theory in powers of the fluctuating order parame- 
ter and potential and further averaging over them using 
the Gaussian approximation 8 or (ii) using the diagram- 
matic technique. For our purpose we choose the diagram- 
matic approach. All diagrams (before impurity averag- 
ing) which contribute to the suppression of the transition 
temperature in Eq. [7| are shown in Fig. 1 . One can see 
that there exist two qualitatively different classes of di- 
agrams. First, the diagrams a) - c) describe corrections 
to the transition temperature due to Coulomb repulsion 
and represent the so called fermionic mechanism of the 
suppression of superconductivity. The second type, dia- 
gram d) , describes a correction to the transition temper- 
ature due to superconducting fluctuations and represents 



a) 



c) 



j — rz 



b) 



d) 



FIG. 1: Diagrams a) - c) describe the correction to the super- 
conducting transition temperature due to Coulomb repulsion. 
The diagram d) describes correction to the transition temper- 
ature due to superconducting fluctuations. All diagrams are 
shown before averaging over the impurities. The solid lines 
denote the electron propagators, the dashed lines describe 
screened Coulomb interaction and the wavy lines describe the 
propagator of superconducting fluctuations. 



the bosonic mechanism. It may seem surprising that we 
classify the diagram (c) as belonging to the fermionic 
mechanism, since this diagram contains both Coulomb 
and Cooper pair propagators. The reason is that, as we 
will show below (see also Ref. l2(jlh there are dramatic 
cancellations between contributions of diagrams of the 
types (a,b) and (c). It is this cancellation that is respon- 
sible for the smallness of the contribution of the fermionic 
mechanism at high temperatures T > T. The diagrams 
of type (c) were not taken into account in Ref. l2ll where a 
different result for the suppression of the transition tem- 
perature was obtained^. In what follows we consider 
both mechanisms of the suppression of superconductiv- 
ity in details. 



A. Suppression of superconductivity due to 
fluctuations of the order parameter: bosonic 
mechanism 



In this section we consider the suppression of the su- 
perconducting transition temperature in granular met- 
als due to fluctuations of the order parameter (bosonic 
mechanism). We will use the diagrammatic technique 
developed in Refs. l5llol The main building block of the 
diagrams that will be considered in this section is the 
Cooperon propagator defined by the diagrams shown in 
Fig. [3 a. In the regime under consideration all charac- 
teristic energies are much less than the Thouless energy 
E T = D/a 2 , where D is the diffusion coefficient. This 
allows us to use the zero dimensional approximation for 
a single grain Cooperon propagator Cq 1 = r|fi n |. The 
resulting expression for the Cooperon is 



C(fi B ,q) 



l (\n n \ + £tl 6y 



(8) 



where q is the quasi-momentum and £l n is the bosonic 
Matsubara frequency. The parameter e q in the right 
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hand side of Eq. JSJ appears due to the electron tun- 
neling from grain to grain, it was defined in Eq. (|ld|) . 

The propagator of superconducting fluctuations, 
K(fl n q) is defined by the diagrams shown in Fig. [21 b 
and c. They result in the following expression 



(c) 



FIG. 2: Diagrams (a) define the Cooperon propagator (shaded 
rectangle), Eq. ©, in terms of the single grain Cooperon 
Co- Diagrams (b) describe the renormalization of the BCS 
interaction vertex due to impurities. The superconducting 
propagator, Eq. @, is represented by the thick wavy line 
and is defined by the diagrams (c) where the thin wavy line 
denotes the bare superconducting propagator. The solid lines 
denote the propagator of electrons and the tunneling vertices 
are denoted as the circles. 



b) 



FIG. 3: Diagrams describing correction to the transition tem- 
perature due to superconducting fluctuations (bosonic mech- 
anism). The diagrams were obtained after averaging the di- 
agram d) in Fig. 1 over the disorder. The solid lines de- 
note the propagator of electrons, the dotted lines describe the 
elastic interaction of electrons with impurities and the wavy 
lines describe the propagator of superconducting fluctuations. 
The shaded rectangle and triangle denote the Cooperon, see 
Eq. JHJ, and impurity vertex of granular metals respectively. 
The tunneling vertices are denoted as the circles. 



#(fi«,q) 



T r v V 2 4ttT 



with tp( x ) being the digamma function. 



(9) 



The diagrams describing the correction to the tran- 
sition temperature in the lowest order with respect to 
the superconducting fluctuation propagator K(S} n , q) are 
shown in Fig. [3] 

Deriving the analytical result for the diagrams in Fig. [3] 
it is important to take into account the fact that the sin- 
gle electron propagator itself gets renormalized due to 
electron hopping. Tunneling processes give rise to an ad- 
ditional term in the self-energy part of the single electron 
propagator, see Fig. 0] 



2dg T 6, 



(10) 



where To is the unrenormalized electron mean free time. 
Although the second term in the right hand side of 
Eq. QlOfl is much smaller than the first one, it is im- 
portant to keep it because the leading order contribution 
in Tq 1 to the correction to superconducting transition 
temperature cancels. 

The contribution of each diagram in Fig. [21to the sup- 
pression of superconducting critical temperature is pre- 
sented in Appendix A. Here we write the final expression 
for the contribution of the diagrams (a)-(f) in Fig. 01 to 
the superconducting transition temperature 



l>0 



K(Q n ,<j)[2\e n \ + \2s n - Sl n \ + e q £] 
4(|2e„ -n n \ +e q (5) 2 



E 

»(e«-n„ 



l<0 



£ 2(|0„| + £ q 5) 



ill) 



Here the summation is going over the quasi-momcntum, 
q and over the fermionic, e n = TrT(2n + 1) and bosonic, 
Q n = 2irTn Matsubara frequencies. 

The main contribution to the suppression of supercon- 



ducting transition temperature in Eq. (|ll|l comes from 
the region of classical fluctuations and is given by the 
term with Q n = 0. Performing summation over the 
fermionic Matsubara frequency e n in Eq. (|11|) we obtain 
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the following result 



/AT C 



14C(3) 



Eh 



(12) 



where is the zeta-function. Equation (|12H for the 
suppression of superconducting transition temperature is 
valid outside the Ginzburg region otherwise the lowest 
order approximation in the superconducting propagator 
which we used to derive Eq. (|12fl is not justified. Perform- 
ing summation over the quasi-momentum q in Eq. i|12|) 
we obtain the final result for the suppression of super- 
conductivity due to the bosonic mechanism, Eq. Ijlel ). 
The singularity in two dimensional case should be cut at 



momenta qf, 



(Tc/g^S) in accordance with the 



expression for the Ginzburg number 



Gi 



1 

1 Te_ 



2d, 
3d. 



(13) 



The divergence of the correction to the transition temper- 
ature in 2d, Eq. (|12fl . means that fluctuations destroy 
the superconducting long range order, which is to be re- 
covered by introducing the artificial cutoff q m i n . Then 
the critical temperature which we calculate should be 
viewed as a crossover temperature rather than the tem- 
perature of a true phase transition. However, since ex- 
perimentally such a temperature marks a sharp decay of 
the resistivity, the notion of the transition temperature 
still makes a perfect sense. 

It follows from Eq. I|13[) that for 3c? granular metals the 
Ginzburg number is small in comparison with the right 
hand side of Eq. (|12fl and the Gaussian approximation 
for bosonic mechanism is justified for temperatures T c < 
g%,5. In 2c? case the result (11211 holds with the logarithmic 
accuracy in the same temperature interval g^,S >> T c . 

The correction to the transition temperature, Eq. (|12fl . 
can be interpreted as a contribution of the fluctuations 
of the superconducting order parameter. These fluctua- 
tions can be considered as a virtual creation of Cooper 
pairs which are bosons. That is why we call this mecha- 
nism of the suppression of the superconductivity bosonic. 
The correction AT C , can be rather easily obtained via the 
Ginzburg-Landau expansion in the order parameter A (r) 
near the critical temperature. For the granular system 
the free energy functional F [A] can be written as follows 



T i . 2 , . ,4 

^AdV-IA/ 



-^TJ^-A/, (14a) 



where the coefficients f and b are given by 
T-T c 7C(3) 



(14b) 



b) 



FIG. 4: Diagrams describing self-energy corrections to the 
single electron propagator due to a) elastic interaction of elec- 
trons with impurities and b) electron hopping. The solid lines 
denote the bare propagator of electrons and the doted line 
describes the elastic interaction of electrons with impurities. 
The tunneling vertices are described by the circles. 



Eq. (|14afl we obtain the propagator of (A q A_ q ) Q in the 
form 



(A q A_ q ) = - 



T,5 



f + (ttS/8T c ) e c 



(15) 



The correction to the transition temperature can be 
found by calculating the first order in b contribution to 
the self energy T,^: 



<) 



f + (7r$/8T )e q 



(16) 



This correction renormalizes the critical temperature. 
Putting f = in Eq. (|16fl we come to the result ex- 
pressed by Eq. ■ 

The fermionic mechanism of the suppression of the su- 
perconductivity is more complicated and we consider it 
in the next section. 



B. Suppression of superconductivity due to 
Coulomb repulsion: fermionic mechanism 

In this section we consider the suppression of the su- 
perconducting transition temperature in granular metals 
due to combine effects of Coulomb interaction and disor- 
der (fermionic mechanism). The Coulomb interaction in 
granular metals is screened by surrounding electrons as 
in any metal. The diagrams that describe the screened 
effect are presented in Fig. EI The diagrams (a) define the 
diffusion propagator in granular metals. As in the case 
with the Cooperon, Eq. (JSJ, we can consider the single 
grain Diffusion in the zero dimensional approximation 
such that Dq 1 = r|O n | and for the diffusion propagator 



we obtain the expression 
D(Q n ,q)=r- 



EqSy 



(17) 



that coincides with the Cooperon propagator in Eq. JSJ. 
The diagram (b) in Fig. [S] describes the renormalization 
of the Coulomb vertex due to impurities and the diagram 
(c) defines the screened Coulomb interaction 



with Ji t i±i = Jq = 7T5 , t/16T' c being the Josephson cou- 
pling between the grains. Neglecting the quartic term in 



V(fi„,q) 



C(q) 



+ 



2e c 



0, 



£ q <5 



(18) 
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FIG. 5: Diagrams (a) define the Diffusion propagator 
D(fl n ,q), Eq. (1171 . in terms of the single grain diffusion 
propagator Do(£ln)- Diagrams (b) define the renormaliza- 
tion of the Coulomb vertex due to impurity scattering. Di- 
agrams (c) define the dynamically screened Coulomb inter- 
action, Eq. 1181 . The thin dashed lines represent the bare 
Coulomb interaction while the thick lines denote the screened 
Coulomb interaction. 

where C(q) is the Fourie transform of the capacitance 
matrix which has the following asymptotic form at q <C 
a" 1 



The correction to the critical temperature due to 
Coulomb interaction before averaging over impurities is 
given by the diagrams a-c in Fig 1. Averaging over the 
impurities leads to rather complicated formulae. We will 
consider the contributions from the diagrams a,b and c 
separately presenting the total correction to the critical 
temperature due to Coulomb interaction as 

=(X 1 ) + (X 2 ), (20) 

where the term X\ represents the contribution of the 
diagrams a and b in Fig. 1 while X 2 represents the con- 
tribution of the diagram c and (...) means averaging over 
the disorder. The two terms in the right hand side of 
Eq. (|2UI) have a transparent physical meaning: the term 
X\ describes the renormalization of the Cooperon due 
to Coulomb repulsion while the term X 2 describes the 
vertex renormalization. After the disorder averaging the 
terms (Xi) and (X 2 ) are represented in Figs. El and [7| 
respectively. The evaluation these diagrams is presented 
in Appendix B. The final expression for the correction to 
the critical temperature due to the fermionic mechanism 
has the following form 



/AT C 



\ T 



-4T J2 V(n n ,<i) 

q,f2„>0 

2F(n n )K(n n ,q) 



F(n n 
M) 2 



£q(5 



(fi„ + e q <5) 2 fi„ 4ttT fi„ - ( £q <5) 2 



2nT J 



[fi„ - ( £q <5) 2 



(V>[ 1/2 + (fi„ + e^/AnT] - fj>[ 1/2 + fi„/27rT]) 



(21a) 



where we introduced the notation 



F(fi„) = V(l/2 + fi„/27rT) - ^(1/2), 



(21b) 



and the propagator K(Cl n ,q) was defined in Eq. 10. 
The summation in Eq. I|21a|) is going over the quasi- 
momentum q and the bosonic Matsubara frequencies fi„. 
The propagator of the screened electron-electron interac- 
tion, V(fi„,q), Eq. (|18|) . in the limit when the charging 
energy, E c , is much larger than the average mean level 
spacing, S, can we written as 



V(fi»,q) = 



2£ c (q)(|fi„|+£ q 5) 



4s q ^c(q) 



in. 



(21c) 



where £c(q) = e 2 /2C(q) is the charging energy. 

An important feature of Eq. Ij21a(l is that the expres- 
sion in the big square brackets in the right hand side 
vanishes at q — > for any frequency fi„. This is due 
to the fact that the potential T^(fi„,q = 0) represents 
a pure gauge and thus it should not contribute to the 
thermodynamic quantities such as the superconducting 
critical temperature. As a consequence one can see that 



the contribution of the frequencies fi„ that belong to the 
interval 



e q <5 < fi„ < e q £:c(q) 



(22) 



have a small contribution to the critical temperature cor- 
rection given by Eq. (|21a|l . 20 . For this reason we do not 
expect the Coulomb energy Eq to appear in the final re- 
sult. At the same time the frequencies fi„ that belong 
to the interval l|22l) are fully responsible for the loga- 
rithmic renormalization of the integranular conductance 
g T (T) = g T - (I /2nd) \n(g T E c /T), where the Coulomb 
energy explicitly appears in the result. 

The expression in the r.h.s. of Eq. I|21a|l is quite com- 
plicated; we cannot derive a simple result at any arbitrary 
temperature. Further on we will consider only the lim- 
iting cases T > r and T < T where the calculations are 
considerably simplified. 

If the temperature T is sufficiently small, T -C T, the 
summation over the Matsubara frequencies can be re- 
placed with the integral. One can easily see that the 
singularities at fi = e q <5 in the second and third terms in 
Eq. (|21a|l cancel each other. Their appearance, in fact, is 
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A A 



d) 



A 



g) 



T 
i 



b) 



~i 



FIG. 6: Diagrams obtained from diagrams a and b in Fig. 1 
after disorder averaging. The solid lines denote the propaga- 
tor of electrons, the dashed lines describe screened Coulomb 
interaction and the dashed-dot lines describe the elastic inter- 
action of electrons with impurities. The shaded rectangle and 
triangle denote the renormalized Cooperon, see Eq. JSJ, and 
impurity vertex of granular metals respectively. The indices 
i and j stand for the grain numbers. The tunneling vertices 
are described by the circles. 



from the low momenta q <C 1/a where the energy e q 
can be written as e q = <7rq a , and the granular system 
becomes equivalent to a homogenously disordered one. 
Summation over q with the logarithmic accuracy leads 
to the final result lflc|) for 2D case. No wonder that 
the result, Eq. (fTcfl in 2 D agree s with the known result 
for disorder metals Refs. l8l9lllljl2CI Equation (|Tc|> in the 
2D case has a universal form and is expressed in terms 
of the tunneling conductance g T and the effective relax- 
ation time (g T (5) _1 . For the homogenously disordered 
samples the latter time should be replaced by the elastic 
scattering time r. 

In the opposite limit, T ^> T the quantity e^S can 
be neglected with respect to the Matsubara frequencies 
f2„ and the result is drastically different from the one 
given by Eq. (|Tc|) . The potential V(f2 n ,q) in this limit 
takes the form V(f2„,q) = |£l„|/2£ q such that e q cancels 
in the main approximation. Summation over Matsubara 
frequencies then leads to the correction Eq. Ijlbfl . 

One thus can see from the above that in the limit 
T > T the fermionic mechanism of the suppression of 
the superconductivity is no longer efficient. This can be 
seen rather easily in another way using the phase ap- 
proach of Ref. 6. Following these works one decouples 
the Coulomb interaction, Eq. I|5c|) by integration over a 
phase <j) 



a pure artifact of the representation of the result in terms 
of the ip— functions. With logarithmic accuracy one can 
leave only the first term in Eq. (|21a[l : this results in 



a) 



(f 



-E 1 

2tt ^ e Q 



ln 2 ^, 
T 



e q 5 > T. 



(23) 



In the 3d case one can neglect the q— dependence un- 
der the logarithm in Eq. (23), and the summation over 
quasimomentum leads to the logarithmically accurate fi- 
nal result (|lcf) . In two dimensions, the main contribution 
in summation over quasimomentum in Eq. 1)231) comes 



w = 



A 



b) 



FIG. 7: Diagrams describing vertex renormalization obtained 
from the diagram c) in Fig. 1 after averaging over the disorder. 
All notations are the same as in Figs. 2 and 4. 



exp 



e 
~2 



J exp I -i J n i( T )^>i ( T ) dT - Yl J 



'/^(t)^-(t) j Do. (24) 



where n^r) = ftp* (ri, r) ip (r^, r) dr, is the electron 
number in the i-th grain and ip are fermionic fields. With 
this decoupling the variable <j) plays a role of an additional 
chemical potential. In the limit T ^>> T (and only in this 
limit) the phase <j> can be gauged out via the replacement 



ip (n, t) ->• ip (t u t) exp [-ifc (t)] 



(25) 



Substituting Eq. l(2*K|l into Eq. (|5al - I5djl (or to be more 
precise, into the corresponding Lagrangian in the func- 
tional integral representation) we immediately see that 
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the phase <f> enters the tunneling term, Eq. I|5d|l only. 
However, this term is not important in the limit T>T 
and we conclude that the long range part of the Coulomb 
interaction leading to charging of the grains is completely 
removed. Therefore, the effect of the Coulomb interac- 
tion on the superconducting transition temperature must 
be small and this is seen from Eq. (|lb|) . This conclusion 
matches well the fact that the upper limit in the loga- 
rithms in Eq. ijlc|) is just Y and at temperatures exceed- 
ing this energy the logarithms should disappear. At low 
temperatures T < Y the phase description does not apply 
and the non-trivial result, Eq. I|lc[) . appears. This result 
is of the pure quantum origin, and interference effects are 
very important. In the limit of high temperatures T > Y 
the interference effects are suppressed, that is why the 
fermionic mechanism of the suppression of the supercon- 
ductivity is no longer efficient. 



VI. DISCUSSION 

We have described the suppression of the supercon- 
ducting transition temperature due to (i) fluctuations 
of the order parameter (bosonic mechanism) and (ii) 
Coulomb repulsion (fermionic mechanism) in granular 
metallic systems at large tunneling conductance between 
the grains, g T ^S> 1. We have calculated the correction 
to the transition temperature for 3c? granular samples 
and films. We have demonstrated that at temperature 
T > g T S the suppression of superconductivity in granu- 
lar metals is determined by the bosonic mechanism while 
at low temperatures, T < g T S, the suppression of super- 
conductivity is dominated by the fermionic mechanism. 

The bosonic mechanism has a classical origin. The 
fermionic mechanism is of the quantum origin and is 
relevant at low temperature T < Y where quantum in- 
terference effects are pronounced. In the opposite limit 
T > Y the coherence is lost and the fermionic mechanism 
of suppression of the conductivity is no longer efficient. 
Thus the classical bosonic mechanism is more important 
at high temperatures, whereas the quantum fermionic 
one is efficient at low temperatures. 

The results of the theoretical study presented in this 
paper can be checked experimentally. Apparently, the 
best way to do this is to study the dependence of the 
superconducting transition temperature T c of granular 
metals as a function of the tunneling conductance g T . 
Practically this can be done by studying several granular 
samples with different tunneling conductances (different 
oxidation coating). The experimental curves for T c sup- 
pression should have a different slope at high T c > g T S 
and low T c < g T 8 critical temperatures due to the fact 
that the suppression of the superconductivity is given by 
the two different mechanisms. Moreover, the difference 
between the dependence of the critical temperature on 
the tunneling conductance T c (g T ) in granular and homo- 
geneously disordered systems can give an information on 
the morphology of the sample, i.e. answer the question 



whether the samples are homogeneously disordered or 
granular. 

Another interesting consequence of our results in 
Eqs. (JIJ is the following: since at low critical temper- 
atures T c < g T S the suppression of superconductivity 
in granular metals is given by the fermionic mechanism 
and upon the substitution g T 5 — > t -1 it coincides with 
the proper result for homogeneously disordered sam- 
ples^ one can generalize the renormalization group re- 
sult by FinkelsteiniS for the T c suppression. The latter 
result obtained for homogeneously disordered films can 
be applied to the case of the granular superconductors 
upon the proper substitution for the diffusion coefficient 
D = g T Sa 2 , where a is the size of a single grain. 
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APPENDIX A: EVALUATION OF DIAGRAMS 
FOR THE BOSONIC MECHANISM 

In this appendix we evaluate diagrams for bosonic 
mechanism presented in Fig.[3]beginning with considera- 
tion of the diagrams a), b) and c). These three diagrams 
can be conveniently combined in a single diagram intro- 
ducing the Hikami bo3*2£ as it is shown in Fig. [S] For zero 
dimensional grain (all characteristic energies are much 
less than the Thouless energy) Hikami box is given by 
the following expression 

H(e ni ,e n2 ,e n3 ,s ni ) = 2 7 rr 2 (|e„ 1 | + \e n2 \ + \e n3 \ + \e ni \), 

where e ni are fermionic Matsubara frequencies. Using 
Eq. (|A1|) and evaluating the diagram shown in Fig. 6b 
we obtain the following result for the sum of these three 
diagrams 

^ T 2c \ " 3|e w | + \fl n - gnj T/, n \ 

e n (s n -n n )>0 n[l " 

(A2) 

where K(Q n ,q) is the propagator of superconducting 
fluctuations defined in Eq. ©. 

The sum of the diagrams d) and e) in Fig. |3| is given 

by 

The above expression for diagrams d) and e) in Fig.[3Im 
fact includes an additional contribution coming from the 
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k>3 



C1J4 





b) 




FIG. 8: Diagrams a) describe "bosonic" Hikami box, 
Eq. lAH . Using Hikami box the sum of three diagrams a), b) 
and c) in Fig. [3] can be conveniently represented as a single 
diagram shown in Fig. 6b. All notations are the same as in 
Fig. 3. 



diagrams a), b) and c) in Fig. [3] that was not included 
in Eq. I|A2J| . This additional contribution appears due 
to the fact that the single electron Green function self- 
energy has a correction resulting from the renormaliza- 
tion of the Green function self-energy due to intergran- 
ular tunneling, see Eq. 1)10(1 and Fig. 4. This self energy 
correction is negligible in the diffusive limit (r — ► 0), nev- 
ertheless it gives a finite contribution to the sum of the 
diagrams a) - c) because each of these diagrams diverge 
in the diffusive limit as 1/r while their sum remains fi- 
nite due to cancellation of the leading orders in 1/t. The 
additional finite term appears because the diagrams b) 
and c) have impurity lines that are determined by the 
bare mean free time to while in all other places the mean 
free time appears through the Green function self energy 
that contains the renormalized r. This additional con- 
tribution could have been written as an extra constant 
term St^dgT^ in the Hikami box. We, however, find it is 
natural to "redirect" this term to the diagrams d) and 
e) since these diagrams are also proportional to the tun- 
nelling conductance gx- 

Finally, for the diagram f) in Fig. we obtain 

^ K(Q n ,q) 1 



ttT 2 S 



„-n„)<o 



\n r , 



S a 5' 



(A4) 



Adding all the contributions given by Eqs. (|A2|) - (IA4I) 
we obtain the correction to the superconducting tran- 
sition temperature due to bosonic mechanism presented 
in Eq. (JTTJ- 



the right hand side of Eq. 12011 that can be written as a 
sum of two terms 



(X t ) = <A7) + 



(Bla) 



where (X°) and (X^) represent the contributions of the 
diagrams a-f and g-j in Fig. respectively. The sum of 
the diagrams a), b) and c) can be presented as a single 
diagram with the help of the Hikami box shown in Fig.|5^, 
exactly as in the case of the bosonic diagrams considered 
in Appendix A. The corresponding Hikami box shown in 
Fig-IHt 1 differs from the "bosonic" Hikami box only by the 
arrow directions and is given by the same Eq.|JlJ The 
sum of the diagrams a), b) and c) in Fig. HJ] thus results 
in the following expression 



-2ttT 2 



E 



,(£„-n„)<o 



(|e n | + |gnj/2) 
el(\n n \+e^ 



V(Q n ,q), (Bib) 



where summation is going over the quasi-momentum, q, 
fermionic, e„ and bosonic, f2„ Matsubara frequencies. 
The propagator of the screened electron-electron inter- 
action, V(fl n ,q) in Eq. l(Blb|l was defined in Eq. (|21c(l . 

The sum of the diagrams d) and e) in Fig. results in 
the following contribution 

E £ 2(\ n \ +£ 6) 2 £ ^ ( B1C ) 

while the diagram f) is given by 



e„(e„-n„)>0 



e\ (|2£ n -fi„| + £q <5) 



(Bid) 



Summing up all contributions in Eqs. ljBlb|) - (|Bld|l we ar- 
rive to the following expression that represents the con- 
tribution of diagrams a)-f ) in Fig. HJ 






I 



"i 



"T- 



b) 




APPENDIX B: EVALUATION OF DIAGRAMS 
FOR THE FERMIONIC MECHANISM 

In this appendix we evaluate diagrams for fermionic 
mechanism presented in Figs. El and [7| . We begin our 
analysis with the evaluation of the contribution {Xi) in 



FIG. 9: Diagrams a) describe "fermionic" Hikami box. Using 
Hikami box the sum of three diagrams g) , h) and i) in Fig. |S| 
can be conveniently represented as a single diagram b). All 
notations are the same as in Fig. HJ . 
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TO 



2 E 

q 



E 

£„(£„-O„)<0 



2^(O n ,q) 



y(n n ,g) 

kn|(|fin|+e q <5) 2 el(\n n \ + e^S) 



(Blel 



Now we turn to the evaluation of the diagrams g-j in 
Fig-EDthat are represented by the term {X\) in Eq. IjBlajl . 
Calculation of the sum of the diagrams g), h) and i) 
again can be reduced to the evaluation of a single dia- 
gram shown in Fig. 7b resulting in 



-7TT 2 



E 



|fi„| v(n n , q ) 



(\fln\+ £q£) 2 M|£n - On I 



(Blf) 



e„(e„-f!„)<0 

The diagram j) in Fig. El is given by the following expres- 



sion 



-ttT z 



E 



£q 5 v(n n ,q) 



e n (e n -Q n )<0 



(|0„| +e q( 5) 2 |e„||e r - r 



(Big) 



while the diagram k) results in 



E 



V(n n> q) 



(\2e n - fi n | +£q(5)|£ n ||£ n - Cl n \ 

(Blh ) 

Summing up the above contributions, Eqs. (|Blfl) - (|Blh(l . 
we obtain the expression representing the sum of the di- 
agrams g) - k) in Fig. HJ1 



(Xl) = -nT 2 



q 



E 

e„(e„-O„)<0 



n n \(\n n . 



E 



v(n n ,q) 



e„(e„-O„)>0 



£« £r, 



I ( I 



e q <5) 



(Bli) 



Now we turn to the evaluation of the vertex renormal- 



malization of the effective interaction vertex between the 



ization which is given by the term (X 2 ) in the right hand Coulomb and Cooper pair propagators. The renormal- 
side of Eq. (|2"U)l . The corresponding diagrams are shown ized vertex r(f2„) is given by the sum of two diagrams 
in Fig. Averaging over impurities results in the renor- shown in Fig. 03 that lead to the following expression 



r(a,) = 



T 



a. 



e a S 



o<e„<n r , 



£„>0 

r 



■rt °n 



O / If 



£q<5 



(Bij) 



Using Eq. | |Blj| ) the resulting expression for the term over the fermionic frequencies in Eqs. i|BlcllBlijl in terms 
(X2) in Ea. (|20|l can be written as of the di-gamma functions after some rearrangements of 

different terms in Eqs. (|Blel IBlifl we obtain Eq. I|21a|) 
(X 2 ) — 8tt 2 T V" V(n n ,q) K(S} n ,q)T 2 (tt n ), (B2) for the suppression of superconducting transition tem- 



q,f2„>0 



perature due to fermionic mechanism. 



where if(f2 n ,q) is the propagator of superconducting 
fluctuations defined in Eq. Expressing summations 
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